Numerical optimization of longitudinal collimator geometry for novel x-ray field.

Objective. A novel x-ray field produced by an ultrathin conical target is described in the literature. However, the optimal design for an associated collimator remains ambiguous. Current optimization methods using Monte Carlo calculations restrict the efficiency and robustness of the design process. A more generic optimization method that reduces parameter constraints while minimizing computational load is necessary. A numerical method for optimizing the longitudinal collimator hole geometry for a cylindrically-symmetrical x-ray tube is demonstrated and compared to Monte Carlo calculations.Approach. The x-ray phase space was modelled as a four-dimensional histogram differential in photon initial position, final position, and photon energy. The collimator was modeled as a stack of thin washers with varying inner radii. Simulated annealing was employed to optimize this set of inner radii according to various objective functions calculated on the photon flux at a specified plane.Main results. The analytical transport model used for optimization was validated against Monte Carlo calculations using Geant4 via its wrapper, TOPAS. Optimized collimators and the resulting photon flux profiles are presented for three focal spot sizes and five positions of the source. Optimizations were performed with multiple objective functions based on various weightings of precision, intensity, and field flatness metrics. Finally, a select set of these optimized collimators, plus a parallel-hole collimator for comparison, were modeled in TOPAS. The evolution of the radiation field profiles are presented for various positions of the source for each collimator.Significance. This novel optimization strategy proved consistent and robust across the range of x-ray tube settings regardless of the optimization starting point. Common collimator geometries were re-derived using this algorithm while simultaneously optimizing geometry-specific parameters. The advantages of this strategy over iterative Monte Carlo-based techniques, including computational efficiency, radiation source-specificity, and solution flexibility, make it a desirable optimization method for complex irradiation geometries.

Technical note: Dose voxel resolution determination for Monte Carlo electron beam simulation in thin targets.

BACKGROUND: Monte Carlo particle simulation has become the primary tool for designing low-energy miniature x-ray tubes due to the difficulties of physically prototyping these devices and characterizing their radiation fields. Accurate simulation of electronic interactions within their targets is necessary for modeling both photon production and heat transfer. Voxel-averaging can conceal hot spots in the target heat deposition profile that can threaten the integrity of the tube. PURPOSE: This research seeks a computationally-efficient method of estimating voxel-averaging error in energy deposition simulations of electron beams penetrating thin targets to inform the appropriate scoring resolution for a desired accuracy level. METHODS: An analytical model to estimate voxel-averaging along the target depth was developed and compared to results from Geant4 via its wrapper, TOPAS. A 200 keV planar electron beam was simulated to impinge tungsten targets of thicknesses between 1.5- and 12.5- µ m ${{\umu {\rm m}}}$ . For each target, the model was used to calculate the energy deposition ratio between voxels of varying sizes centered on the longitudinal midpoint of the target. Model-calculated ratios were compared to simulation outputs to gauge the model's accuracy. Then, the model was used to approximate the error between the point value of electron energy deposition and a voxel-based measurement. RESULTS: The model underestimates error to within 5% for targets less than 7.5- µ m ${{\umu {\rm m}}}$ in thickness with increasing error for greater thicknesses. For the 1.5- µ m ${{\umu {\rm m}}}$ target, calculations of the point-vs.-voxel energy deposition show an 11% averaging effect between the midpoint and a 1.5- µ m ${{\umu {\rm m}}}$ voxel. Energy deposition profiles along the target depth were also calculated in the Monte Carlo for reference. CONCLUSION: A simple analytical model was developed with reasonable accuracy to guide Monte Carlo users in estimating the appropriate depth-voxel size for thin-target x-ray tube simulations. This methodology can be adapted for other radiological contexts to increase robustness in point-value estimations.

Comparison of novel shielded nasopharynx applicator designs for intracavitary brachytherapy.

PURPOSE: Nasopharyngeal brachytherapy is limited in part by the radiotolerance of nearby organs like the soft palate. This study explores several novel shielding designs for an intracavitary applicator to significantly reduce soft palate dose while adhering to the constraints of standard treatment procedure. METHODS: The Monte Carlo code TOPAS is used to characterize each prototype under typical high-dose-rate treatment conditions. Mucosal surface dose maps are collected to evaluate the shields on their dose reduction to the central and soft palate planning points and uniformity in their shielding profile. Practicality with respect to patient comfort and pretreatment imaging is discussed. History-by-history standard deviations are calculated for each simulation. RESULTS: A design with elliptical tubing containing bundles of tantalum wires provides the most significant attenuation with 39% and 27% dose reduction to the center and soft palate locations, respectively. Another design utilizing miniature lead spheres loaded into a constructed cavity shows 27% and 24% dose reduction to the same locations while providing more uniform shielding and several practical benefits. Both shields are designed to be completely removable for applicator insertion and pretreatment imaging. The mean and maximum standard error of relative dose measurements was 0.36 and 1.14 percentage points, respectively. CONCLUSION: Each shielding design presented in this study provides a novel approach to safely and effectively shield healthy tissue during intracavitary nasopharyngeal brachytherapy. Analysis performed using Monte Carlo suggests that the design using metal spheres most practically shields the soft palate and should be advanced to the next stages of clinical optimization.

Paradigm Shift in Radiation Treatment Planning Over Multiple Treatment Modalities.

The inverse planning simulated annealing optimization engine was used to develop a new method of incorporating biological parameters into radiation treatment planning. This method integrates optimization of a radiation schedule over multiple types of delivery methods into a single algorithm. We demonstrate a general procedure of incorporating a functional biological dose model into the calculation of physical dose prescriptions. This paradigm differs from current practice in that it combines biology-informed dose constraints with a physical dose optimizer allowing for the comparison of treatment plans across multiple different radiation types and fractionation schemes.